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Abstract 

We note that the phenomenon of perturbative saturation leads to transverse 
momentum broadening in the spectrum of partons produced in hadronic collisions. 
This broadening has a simple interpretation as parton level Cronin effect for systems 
in which saturation is generated by the "tree level" Glauber-Mueller mechanism. 
For systems where the broadening results form the nonlinear QCD evolution to high 
energy, the presence or absence of Cronin effect depends crucially on the quantitative 
behavior of the gluon distribution functions at transverse momenta kt outside the 
so called scaling window. We discuss the relation of this phenomenon to the recent 
analysis by Kharzeev-Levin-McLerran of the momentum and centrality dependence 
of particle production in nucleus-nucleus collisions at RHIC. 
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1 Introduction 



The phenomenon of perturbative saturation has been the focus of intensive study in recent 
years. Since the appearance of the first RHIC data, saturation based ideas have motivated 
several attempts at understanding bulk properties of ultra-relativistic heavy ion collisions 
such as the multiplicity, rapidity distribution and centrality dependence of particle pro- 
duction P3I2]- I n particular the recent work [2] suggests that the saturating properties of 
the nuclear gluon distribution may be responsible for the N pait scaling of charged particle 
multiplicity at 2 < p t < 8 GeV with centrality in RHIC data. It is undoubtedly true that 
saturation effects suppress the number of gluons below the saturation momentum Q s in 
the nuclear wave function. It is then very plausible that the number of produced particles 
at p t < Qs is also suppressed relative to simple perturbative prediction. The value of Q 2 S 
is estimated to be of order 2 GeV 2 for most central collisions at RHIC [I] . The suggestion 
of [2] goes beyond this simple statement and implies that the suppression persists also at 
higher momenta, namely in the so called scaling window Q s < p t < Q^/Qo- 

On the other hand one expects a competing effect of transverse momentum broadening 
due to multiple rescatterings in the final state, the so called Cronin effect. This should 
enhance the number of particles produced in the intermediate momentum range and thus 
works against the saturation argument. The purpose of this note is to show that there is no 
outright contradiction between the appearance of the Cronin effect and the expectations 
based on the saturation scenario and that, in fact, the Cronin enhancement is inherent in 
some realizations of saturation physics. We show here that saturation models based on 
multiple rescattering lead to the relation 

rfiV sat ^ dN pert 

d 2 p t d?p t 

for all Pt > Q s - This relation leads to Cronin enhancement when comparing gluon pro- 
duction in central and peripheral collisions. 

For saturation models based not on multiple rescattering but rather on coherent sup- 
pression at low x the situation is more complicated. Here it is also true that the number 
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of gluons produced in the intermediate momentum range is greater than the prediction 
of the leading order perturbation theory. This overall enhancement is the result of the 
increase of the number of gluons in the nuclear wave function due to low x evolution. The 
transverse momentum broadening is also present. This is manifested by the anomalous 
dimension generated by the low x evolution, so that the decrease of the evolved distri- 
bution with momentum is very slow. However, whether this broadening results in the 
Cronin enhancement is determined by the behavior of saturated gluon distribution at 
high transverse momenta. The momenta which are important are those above the scaling 
window inside which the value of the anomalous dimension is analytically understood. If 
the distribution above the scaling window approaches quickly the leading order perturba- 
tive one, the Cronin effect is indeed generated. However if the distribution continues to 
vary much slower than 1/pt, the (properly normalized) multiplicity of produced gluons is 
always smaller for central collisions than for the peripheral ones. 

In this light, the results of j2] should be understood in the following way. In a 
particular saturation ansatz considered in [2] for p t in the scaling window, ^pr^- as a 
function of centrality is indeed proportional to the number of participants N part . However 
its value is always greater than that for the leading order perturbative prediction for 
the same momentum, even though the latter is proportional to the number of collisions 

4/3 

iV co ii oc N pait . This holds as long as the transverse momentum in question is greater than 
the saturation momentum for the most central collisions considered. Thus according to 
the saturation scenario, the number of produced particles in the intermediate transverse 
momentum range is enhanced and not suppressed relative to the leading order perturbative 
one. This enhancement in the overall production may or may not result in the Cronin 
enhancement when comparing the production rates for central versus peripheral collisions, 
as the evolution is equally effective for all impact parameters. 

Before we proceed further, we wish to make clear the following points. First, to 
calculate the gluon production we use the k t factorized formalism as in ^Ej- Although it 
has not been proved to hold for this process and is likely not to be strictly valid, one may 
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hope as in |T] that it gives a qualitatively reasonable description of gluon production. 
Second, in this note we only address the multiplicity of produced gluons. For comparison 
with experimental data, this quantity has to be convoluted with gluon fragmentation 
functions to convert it to the number of produced hadrons. This introduces additional 
uncertainties related to our limited knowledge of the dynamics of the system between the 
time of production and the time of hadronization [Sj, and the possible medium- dependence 
of parton fragmentation Thus our results have to be paralleled to those of [2] prior 
to convolution with fragmentation functions. 

The simple fc r factorized formula for the gluon yield at central rapidity in a collision 
of identical nuclei [7j used in [2] is 

7-! da A7r 2 a s S AA {b)N c 1 t 2 

E d^b = jvg-i g J d h My ' kt) My ' Pt ~ kt) ■ (2) 

Here S AA (b) is the overlap area in the transverse plane between the nuclei at fixed impact 
parameter b, y is the rapidity difference between the central rapidity and the fragmenta- 
tion region and 4>a{v, kt) is the intrinsic momentum dependent nuclear gluon distribution 
function, related to the standard gluon distribution by 

d{xG A (x, k 2 t )) 

^ k ^= d 2 ktd% ■ o) 

Eq. (J2J itself is an approximation even within the k t - factorization scheme for two reasons. 
First, the gluon distribution is considered to be effectively impact parameter independent 
and taken at some representative impact parameter inside the overlap region. Second the 
rapidity of both distribution functions is taken to be the same rather than integrating over 
the relative rapidity of the two distributions. This amounts to the assumption that the 
gluons are produced locally in rapidity and keep the same rapidity label as in the parent 
distribution function. The first assumption in principle is easily relaxed by considering 
6-dependent distributions, although it makes the calculation considerably more cumber- 
some. The second assumption is more questionable. One certainly expects "migration" in 
rapidity during the interaction, and in particular the relevant rapidities should depend on 
the produced transverse momentum. This is certainly the case in the collinear factorized 



perturbative formalism, where the parent rapidities are taken at x — kt/ys. However in 
the present case the production is not necessarily through the 2^2 process, and thus the 
rapidities are not fixed in the same way. Keeping these caveats in mind, we now consider 
the implications of Eq. (J2J). 

In a purely perturbative approach, neglecting the effects of saturation the gluon dis- 
tribution function at impact parameter b has the shape 

-pert/, 7 \ <Xs(N2 - 1) „ PpartW 1 m 

Va (h,o) = 7T2 K — o — 71 T a u> 4 

27r 2 2 (kt + AqcdY 

where p pa . r t(b) is the density of participants, taken as reference for A — A collisions. The 
role of Aqcd i n the denominator is to regulate the gluon distribution in the infrared. The 
additional numerical multiplicative factor K reflects the fact that the gluons at low x 
originate not only from the valence quarks, but also from the sea quarks and energetic 
gluons. 

2 Models for the gluon distribution 

Saturation effects modify the gluon distribution so that it is suppressed at low momenta 
\kt\ < Q s - There are several models in the literature which provide such saturated distri- 
bution function. In this note we consider two types of models. 

2.1 McLerran— Venugopalan gluon [5J 

The McLerran- Venugopalan model jH] achieves saturation by taking into account the 
Glauber-Mueller multiple scattering effects. The intrinsic glue distribution in this model 
was calculated in [TDJ 

We will take the saturation momentum to be x-dependent following 

«(*^ = ^,G(M/x 2 )^, (6) 
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with G(x,k 2 = 1/x 2 ) being the nucleon gluon distribution. We will take the gluon 
distribution in the nucleon to be of the simple perturbative form [TT] 

xG{x, 1/x 2 ) = K as{N ; ~ 1} In + a) . (7) 

Z7V \ X QCD / 

The small regulator a = -^A — ensures that the saturation momentum stays positive 

X c A QCD 

for x 2 ^> x 2 . For the numerical evaluations we choose x c = 3 GeV -1 , such that for 
momenta k t > 0(1 GeV) the sensitivity on the infrared cut off a is negligible. We take 
Aqcd = 0.2 GeV, and a s = 0.5. The saturation scale Q 2 (b) is obtained from the solution 
of the implicit equation Eq. (JHJ), when evaluated at x 2 = 1/Q 2 {b). At b = we fix 
Q 2 = 2 GeV 2 throughout this paper; this corresponds to K = 1.8 in Eq. ((7J). 

The MV distribution appears in the light cone gauge calculation as the average of 
the gluon number density operator in the state with random distribution of color charges 
distributed with the nuclear density ^U]- On the other hand as shown in in the 
covariant gauge calculation of DIS-like processes it accounts precisely for the rescattering 
of the produced gluon inside the nucleus. The effect of these multiple rescatterings is to 
broaden the gluon transverse momentum spectrum by the amount A/c 2 ~ Q 2 S . Accord- 
ingly, the low momentum part of is suppressed relative to the perturbative gluon, 
the region around k t ~ Q s is enhanced, while at large momenta k t ^> Q s there is no 
appreciable change relative to the perturbative gluon. FigHt shows the ratio of the MV 
gluon to the perturbative one as a function of momentum for fixed coupling constant 
a s = 0.5. The multiple rescatterings do not change the total number of gluons but only 
redistribute the gluon momentum. Thus the gluon distribution G(Q 2 ) calculated with 
is the same as the perturbative one for Q 2 ^> Q 2 S . 

In a DIS-like process with a probe directly coupled to gluons, the intrinsic momentum 
distribution would be directly proportional to the spectrum of final state gluons. The 
Cronin effect is therefore present in the MV gluon ab initio. This has been noticed in 
[12] • The question we are interested in, is to what extent this effect shows up in nuclear 
collisions. 
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Figure 1: a). The unintegrated gluon distribution function, normalized to the correspond- 
ing perturbative one, as function of k t for fixed p part = 3.1 fm~ 2 and for Q 2 S ~ 2 GeV 2 ; 
solid curve: gluon distribution in the McLerran-Venugopalan model jH], dot-dashed curve: 
model (fTU|) - (fT2*j) for the evolved gluon distribution. The dashed curve is for the anoma- 
lous dimension 7 = 0.5. b). The ratio /£P cn P h in Eq. (fTTj) for the gluon distributions 
Eq. ffTUj) (dot-dashed curve) and Eq. (fT3*j) (solid curve). 

2.2 Evolved gluons 

The MV gluon distribution does not contain any evolution in x. One way to introduce 
the x dependence is to adopt the Golec-Biernat-Wusthoff procedure [12], whereby the 
saturation momentum is taken to be energy dependent with the factor x~ A , with A = 
0.2 — 0.3. In this paper we are not going to explore the energy dependence of the spectrum, 
and thus for our purposes the energy independent MV ansatz is sufficient. 

Another type of saturated gluon distribution has been used in The energy de- 
pendence in the RHIC energy range is not large enough to allow one to explore the 
perturbatively predicted x dependence of the distribution function. One can nevertheless 
consider (at both 130 and 200 GeV) as being evolved by the perturbative evolution from 
lower energies. Such an energy evolution leaves a distinctive imprint on the k t dependence 
of the gluon. Although the solution of the nonlinear QCD evolution equation ^3] has 
not been analyzed in great detail, its qualitative features have been discussed in [T^]. It 
has been argued in (T^j that in the wide region of momenta Q s (x) < \k t \ < the 
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evolved distribution behaves as 

~Ql 



with the anomalous dimension 

7 = 0.64. (9) 

A slightly different analysis of [2j based on doubly logarithmic approximation suggests 
7 = 0.5. In either case due to the large anomalous dimension, 0^ LE is very significantly 
enhanced over the perturbative 0^ er in the wide range of momenta. At asymptotically 
large momenta 0^ LE again reduces to the perturbative expression. 

As we will see in the following it is important to know how the distribution behaves 
outside the scaling window. The behavior outside the scaling window is not known an- 
alytically. One expects that at asymptotically large momenta the behavior of 0^ LE is 
perturbative, and therefore 7|fc t |^oo ~~ ► 1- The crossover from the scaling with the anoma- 
lous dimension to the perturbative one can in principle be either sharp at the edge of the 
window, or can be very gradual and slow. To explore the possible differences between the 
fast and slow crossovers we will use two parametrizations of the distribution function. 

To model the function with the fast crossover we take for illustration 

where 

Ql = Ql(b) = K2na*N c ?^, (11) 

and 

l + OMwjkt) f . ( 10 Q 2 s (b) y 

1 +w(k t ) \ (k t + AqcdJV 

Here, Q 2 s (b) = 4Q 2 s (b), consistent with Eq. © at p part = 3.1 fm" 2 . 

The parametrization of <pj^ E (kt) is chosen to be consistent with the required behavior 
at large kt as well as at k t = 0(Q S ) ^Hj- The width of the scaling window in our 
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parametrization is ~ 3Q S . For k t < Q s , the gluon 0^ LE is supposed to saturate or grow 
at most logarithmically. In our ansatz this saturation is ensured by the presence of the 
term Q 2 S in the denominator of Eq. (jlUj) . 

To model the possibility of the slow crossover we will take simply the function with 
fixed anomalous dimension: 

0.64 



1NLES/1, \ _ N c 1 ( Ql 



(**) = -dr-^r 77^7 , (13) 



M*) = 7^3£&. (15) 



Although this function never approaches the perturbative asymptotics, for the purposes 
of numerical evaluation it is indistinguishable from a function with slowly varying ^(k). 

We note that Eq. (|10|) parametrizes the gluon distribution directly in momentum 
space. We found that this is the simplest way to generate an acceptable to be used in 
the framework of Eq. (j2j). Alternatively one could try to define via the frequently used 
relation involving the dipole scattering cross section iV(x) [To] : 

^™=zkk J < i4 > 

where 

(x 2 Qg(x 2 )/4)~ 

1 + (x2Q2( X 2)/4)' 

with 7(x) approaching 7 = 0.64 at small values of 1/Q S > x > Aq CD /Q 2 in a way 
similar to Eq. (|12j) . Although strictly speaking <pA{k t ) is defined by the relation Eq. (|14j) 
only at large k t , naively one could expect it to be also reasonable at k t ~ Q s . However 
it turns out not always to be the case. The Fourier transform in Eq. (JT%j) is n °t only 
sensitive to x ~ l/k t but to all x < l/k t . As a result depending on details of the 
parametrization of 7(2;), we sometimes found "oscillating" unintegrated gluon functions 
which for some momenta were even negative. In general, this indicates that one should 
be very careful using relation Eq. ()14j) . as even for a reasonable N(x) it can produce 
unacceptable </>^ LE (fc t ). 

As seen from Figd the qualitative features of the gluon distribution in the MV model 
(JHJ and the scaling models (fTOj) are to some extent similar. When compared to the 

9 



leading order perturbative distribution Eq. (JH) they show suppression at k t <C Q s and 
enhancement at k t > Q s - In FigQJi, we show the ratio 

R pcrt = Mh)/<f>T{h) ■ (16) 

In the scaling model the enhancement is much more pronounced. It is clear that using 
the same ansatz, but with 7 = 0.5 rather than 7 = 0.64 would make this enhancement 
even greater. 

In relation to the MV gluon the perturbative distribution Eq. (j3J) is the relevant 
distribution to be used to model the peripheral collisions. For the evolved functions a 
more meaningful comparison is with the "peripheral" distribution of the same functional 
form as <pj^ E but with a smaller value of Q s . In Fig^> we plot the ratios 
R periph 0^ LEF(S) (b = 0, k t )/p p3uIt {b = 0) 

NLEF(S) ^L EF (S) (fe = 13fmA)//w(6 = 13fm) W 

with Q s corresponding to p part = 3.1 fm -2 and p part = 0.35 fm -2 , (the values of impact 

parameters 6 = and b = 13 fm) respectively. We observe that the ratios for 0^ LEF and 

0^ LES are very different. While the ratio -R^lef exhibits enhancement for the momenta in 

the scaling window, no such enhancement is seen in R^^. The reason is that with the 

ansatz 0^ LEF as long as kt is outside the scaling window of the peripheral distribution, 

but inside the window of the central one (that is for 3Q s (b = 13fm) < k t < 3Q s (b = 0)) 

the ratio -Rnlef ^ s practically equal to R peTt . This is not the case for 0^ LES because of a 

very slow approach to asymptotic behavior. Instead, for all momenta of interest 

^iph _ ( Q 2 s (b = 0) y 1 _ p piiIt (b=13im) 1 ^ 

^ NLES - \Qi(b= 13 fm) ) ~ pjb = 0) < 1 ' U8j 

We note that the physics of the enhancement of R pcrt for 0^ LE is different from that of 

0^ v . According to jT^j and j2] the anomalous dimension of the evolved distribution 

is a direct consequence of the linear BFKL (or doubly logarithmic) evolution and is 

not significantly affected by the nonlinearity of BK equation. As opposed to multiple 

rescatterings resummed in 0^ v , the BFKL evolution greatly increases the total number 

of gluons relative to the perturbative distribution. Thus the enhancement of 0^ LE is not 
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due to the gluon number conserving redistribution of the gluon momentum, but rather 
due to the BFKL growth of the total number of gluons. The slow fall off of <^{ LE with 
momentum in the scaling window and the associated clear momentum broadening is the 
result of the BFKL diffusion, which fills the phase space very far from the momentum at 
which the distribution is peaked at initial energy. 

Although some qualitative features of saturating gluon distributions appear to be 
model-independent, quantitative results can vary significantly, see Fig^i and Fig^. 
In the rest of this paper, we discuss the implications of this model-dependence for the 
spectrum Eq. (J2J), and we study the behavior of dN/d 2 Pf clS 8j function of transverse 
momentum and centrality. 

3 Gluon production in A-A collisions 

For the perturbative gluon distribution Eq. (0J), the main contribution to the integral in 
Eq. (j2J) comes from the region of phase space k t — Pt ~ and kt ~ 0, pt — k t ~ p t . 
The contribution from the bulk of the phase space k t ~ k t — pt ~ Pt is logarithmically 
suppressed due to the fast decrease of the perturbative distribution function with k t . One 
finds 



dNP eTt (b) 



y -_r 2SMb) w c ^[^ +2 ^' (19) 



dy d 2 p t 

where 7^ is the Euler constant. With Eq. (fTTj) this spectrum scales with N co u oc Saa^) Pp art (fr) 
as expected perturbatively. 

For the saturated gluon distribution in the MV model the gluon yield Eq. (j2J) is 

expressed by 



dN(b) 



dy d 2 p t 



y=0 



Since small momenta pt are suppressed in the MV gluon distribution, expression PT)|) is 

smaller than the perturbative one. For small \p t \, the x-dependence of the saturation scale 

Eq. (jHJ) is frozen, and one finds 

dN[b) Pt «Qs( b ) _ 1 N 2 -lQ 2 s (b)S AA (bl HAy 



dy d 2 p t 
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4?r 3 a s N c p 2 t 



Thus, for small transverse momentum, the spectrum (|21jl scales with the number of par- 
ticipants, i.e. 

N paIt (b) = S AA (b) p paIt (b) . (22) 

Numerically, however, we find that the limit (}2T]) provides a fair approximation for the 
full expression (J2U|l in a very small p part (fe)-dependent region below p t < 0.1 GeV only. 

In FigEl we plot the result of the numerical evaluation of the formula Eq. (|2Tlj) nor- 
malized to the peripheral (perturbative) yield Eq. (|T§|). i.e. 

dN(b) / dN^jb) ^ 
dyd 2 p t I dyd 2 p t 

for Ppart = 3.1 (fm -2 ), and with Q 2 = 2 GeV 2 . This corresponds to the normalized ratio 
of central over peripheral yields 

1 dN(b = 0) / 1 diV(6= 13fm) 

^coii(6 = 0) dyd 2 p t I N coll (b= 13 fm) dyd 2 p t ' ( } 

To understand the p t - dependence of (J23j) qualitatively, we consider the integral / d 2 k t 
4>A{y, k t ) 4>A(y,Pt — k t ) in (J2J). For a saturated gluon distribution and \p t \ > Q s , one does 
not gain a logarithmic enhancement factor from the low momentum region, \k t \ <C \pt\- 
However due to enhancement in the intermediate region one gets a bigger contribution 
to the integral from \k t \ ~ \k t — p t \ ~ Q s . Since both factors of <fi are enhanced in 
this region, we expect that for some momentum range this enhancement will overcome 
the suppression in the small momentum range, and therefore will lead to a net excess 
of produced gluons relative to the perturbative result. This is seen in FigO as a small 
but clear Cronin enhancement of the produced gluon number for momenta just above 
the saturation scale. The amount of enhancement depends on the value of the coupling 
constant, but qualitatively the phenomenon persists for any a s . 

For the evolved distributions </>^ LE ) since the kt dependence is very slow, the contribu- 
tion to the integral comes from a very large range of momenta - for p t in the scaling region 
the integral is dominated by \kt\ ~ \kt~ Pt\ ~ \Pt\- This leads to a significant enhancement 
of gluon production for all \pt\ > Q s relative to the perturbative expression. It is however 
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Figure 2: Cronin effect in the p t - dependence of gluon production yields for head-on A- A 
collisions. Solid curve for the MV-gluon distribution normalized to the perturbative yield. 
The dot-dashed curve is for the evolved gluon distribution (|10[). the dashed one for the 
evolved gluon distribution (|T3|) . both normalized as in Eq. (|24|) . 

more interesting to consider the ratio of the central to peripheral yields. Again, taking 
Ppart(fr = 0) = 3.1 fm~ 2 and p pa , r t(b = 13) = 0.35 fm~ 2 we display the results of the nu- 
merical integration of Eq. (J2J) in Fig|21 This plot mirrors the plot for a single distribution 
FigCJx The distribution </>^ LEF displays a clear Cronin effect similar to the MV gluon, 
while 0^ LES shows uniform suppression for the central/peripheral ratio for all momenta. 
Although the ratio of the yields for (/>^ LEF is smaller than unity at p t ~ 10 — 20GeV, we 
have checked numerically that at very large p t it slowly approaches unity from below. We 
thus conclude that the properties of the ratio are very sensitive to the way in which the 
distribution behaves outside the scaling window. 

We discuss now to what extent the gluon spectrum (j2J) shows an approximate N paxt or 
iVcoii scaling in some kinematic regime. To this end, we plot the gluon yield normalized 
to the yield for peripheral collisions at b = 13.0 fm, corresponding to p part = 0.35 fm~ 2 
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Figure 3: Centrality dependence of gluon production yields (J2J) in A- A collisions as a 
function of p part (fm -2 ), normalized to the yield in peripheral collisions, Eq. (|25[). Curves 
are calculated for the MV gluon distribution (jSJ) and different values of p t : solid curve: 
p t = 0.25 GeV, dashed curve: p t = 1.0 GeV, dot-dashed curve: p t = 3.0 GeV. 

m, 

2 dN(b) J dN{b= 13.0) 
A r P art(&) dyd 2 p t / dyd 2 p t 

This quantity is plotted for the MV gluon in FigJ3]and for the evolved gluon distributions 
(fTUj) and (fT3*j) in FigEJ for various fixed values of p t as a function of p part - We also 
replot in FigElthis quantity for the evolved gluon as a function of N part using the relation 
between p part and A r part given in [T]. We explore the dependence on p part beyond the 
experimentally accessible range to illustrate better the functional dependence. For the 
MV gluon a steep increase with p part indicative of N co a scaling is found for p t > 1 GeV. 
For smaller transverse momentum, e.g. p t = 0.25 GeV, the p part -dependence is seen to 
level off. 

For the evolved gluon the centrality dependence is again sensitive to the large mo- 
mentum behavior. For 0^ LEF at p t = 3 GeV the centrality dependence is similar to the 
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Figure 4: Centrality dependence as in Fig. El but with the gluons of Eq. (fTU|) and Eq. (|13j) 
for p t = 1 and 3 GeV. The solid and dashed lines correspond to Eq. (fTUj) for p t = 1 GeV 
and p t = 3 GeV, respectively. The short-dashed (p t = 1 GeV) and dot-dashed (p t = 3 
GeV) lines are calculated for the gluon distribution (fTB^) . 

MV gluon. One does not recover iVp art scaling even in the enlarged centrality region 
Ppart < 8 fm -2 . The qualitative features of FigsEJ and E] at small and large p paiTt at 
Pt = 3 GeV can be understood as follows: At small p par t, the saturation momentum Q s is 
small so that the scaling window does not exist; the overall yield scales with iV co n, like the 
perturbative one. This iV co ii behavior should persist as long as 3 GeV is above the scaling 
window, 3 GeV > Q^(b)/Qo (for our parametrization Q = 0.5 GeV). This is indeed seen 
for small p part in FigsE] and El In the other extreme, when Q 2 S becomes very large and 3 
GeV lies well inside the scaling window, the main contribution to the yield comes from 
the integration over the momentum in the scaling window and one obtains 

d M«mov- (2 6 ) 

dyd 2 pt 

This dependence is much flatter. For 7 = 0.5 as in j2], it scales with iV part rather than 
N co \\. For the value 7 = 0.64 used here, the dependence is only slightly steeper. This is 
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Figure 5: Normalized yield for evolved gluon distributions as a function of N pait . Legend 
the same as on Fig.4. The largest value iV part = 380 corresponds to p part = 3.1 fm -2 on 
FigH 

the argument given in [2]. However, as we have shown above, at these large values of 
Q s the absolute magnitude of is greater than that of ^^5~ ■ Hence there must 

be an intermediate region of p par t where grows faster than iVcoii so that ^2p t can 

overtake dN d y d 2^ and then stay flat for some region of large p part l . Such behavior is 
indeed seen in FigsEJandEl However, the values of p part for which the curve flattens out, 
are larger than the experimentally relevant ones, see Fig|S| The reason why the scaling 
of Eq. is not reached faster is that the argument leading to Eq. (1261) neglects the 
finite width of the scaling window. Since the function 0^ LEF decreases quite slowly in 
the scaling window, the yield gets significant contributions from momenta above Q 2 s /Qo- 
These momenta bring in additional Q s and therefore p part dependence, and slow down the 
onset of scaling Eq. f)26|). At smaller momenta {p t = lGeV) the normalized yield is flat in 

: Of course g^gs^ does not stay above d ^ d2 ^ for arbitrarily large p par t- When the value of Q s reaches 
p t , dI ^ yd 'ip^ overtakes dy<fip f consistently with our discussion of the previous section. 
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Pp ar t for pp art > 2fm~ 2 , since these gluons are produced below the saturation momentum. 

For 0^ LES on the other hand, the centrality dependence is very flat for all momenta 
we explored, consistent with 

4 Gluon production in p-A or d-A collisions 

Finally, we discuss the gluon production in the situation where the distribution of one of 
the nuclei is perturbative. This situation pertains to proton-nucleus and deuteron-nucleus 
collisions. Following the previous discussion, we use the /c r factorized formula (J2J) for the 
gluon yield, replacing <j>A(y,k t ) 'pAi.ViVt — h) in © by the product of a perturbative 
gluon distribution for the proton and a saturated MV or evolved gluon distribution for 
the nucleus. We also compare the results of this calculation to the gluon production 
cross section derived in the quasi-classical approximation This quasi classical 

expression can be written as 

daPA = 1 ( dhzd 2 y 1 a * CF *- y jr*<x-y) 
d 2 p t dy d?b n J (27r) 2 it x 2 y 2 

(x-y) 2 Q?/4 _ i) + (i _ e- 2 ^/ 4 + 1 - e-^/ 4 )l , (27) 



e 

where b is the impact parameter, which we choose to be b = 0, p t is the gluon's transverse 
momentum, and y its rapidity. Although this expression also has a factorized form in 
momentum space, this form is distinct from Eq. (j2J) [TE] . For the numerical evaluation, 
we regulate the y-integration of the first bracket of (|27j) by an infrared cut-off where 
we choose fx = Aqcd- This cut-off can be introduced using In ^2 = \ J d 2 y ^'(z+y) 2 • 
in Eq. (0), we use for the saturation scale 

Air 2 n N 

Q 2 (x) = — (2Rp) xG(x, 1/x 2 ) , (28) 

c 

where p denotes the nuclear density and xG is the gluon distribution of Eq. (J2J). For the 
numerical analysis, we use 2Rp = p part (6 = 0)/2 = 1.5 fm~ 2 . 

In order to compare with the perturbative behavior, we calculate for central rapidity 
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Figure 6: The gluon yield produced in p-A collisions, normalized to the perturbative 
yield for the MV gluon and to the peripheral yield for the evolved gluon distributions. 
The different curves are for the quasi-classical expression (|27|) normalized to (dashed 
curve), the /c r factorized spectrum with MV gluon (solid curve), and with evolved gluon 
distributions ([10)1 (dot-dashed curve) and Eq. (jT5|) (short-dashed curve), respectively 



y — the asymptotic limit of 1)27)1. 



daP A 



d 2 pt dy d?b 



*a 3 s C F „, nn \ 1 /i r P 2 t 



K (2Rp) - (ln[ 



+ 1lE ~ 1) • 



(29) 



In the limit of small momenta, we quote [3] 



da pA 



d 2 p t dy d?b 



a, Cf 



7T- 



^2 ' 



(30) 



Pt<.Qs " Pt 

which is obtained for "frozen", i.e. constant Q s . We use Eq. (129)1 as the baseline for 
comparison for the MV gluon. For the evolved gluon we calculate the yield using the 
standard factorized formula. We take the proton distribution in the same functional 
shape as that of the nucleus, i.e. Eq. (jT0)l and Eq. (|T5)) with Q s = Aq C d- We compare 
the central yield (p par t = 3.1 fm~ 2 ) to the peripheral yield (p pa rt = 0.35 fm -2 ). Results are 
plotted in FigEl 
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For moderate momenta, i.e. above about twice the saturation scale Q s ~ 1.4 GeV, 
the full rate for the MV gluon shows a Cronin-type enhancement with respect to the 
perturbative one. For small momenta p t < Q s , there is significant suppression, which can 
be immediately deduced when comparing (|30|) with (|29jl. A qualitatively similar behavior 
is found when calculating the gluon spectrum from the factorized ansatz (J2J). As in the 
case of the nucleus-nucleus collisions, we find also that in p-A collisions the Cronin ratio 
in the evolved case depends strongly on the properties of the evolved distribution above 
the scaling window. 

5 Conclusion 

In summary, our study illustrates that quantitative results depend largely on the pre- 
cise model-dependent implementation of saturation effects. A generic qualitative feature 
for the multiple scattering situation is that perturbative saturation leads to Cronin-type 
transverse momentum broadening of the produced gluon spectrum in both nucleus-nucleus 
and proton(deuteron)-nucleus collisions. In the evolved case there is an overall enhance- 
ment of the production yield relative to the perturbative baseline. The spectrum also 
exhibits strong transverse momentum broadening due to relatively large anomalous di- 
mension. On the other hand, the absence or presence of the Cronin effect strongly depends 
on the behavior of the distribution function outside the scaling window. At present the 
shape of the evolved distribution in this momentum range is not known analytically. The 
recent numerical study strongly indicates that crossover from the scaling regime 

to the perturbative one is very slow and gradual, and that the Cronin effect which is 
present in the MV gluon is wiped out by the quantum evolution at high energies. Thus, 
^nles j n (JT3*j) seems to provide a more realistic parametrization of the evolved gluon 
distribution than 0^ LEF in Eq. (jTU|) . 

Although the centrality dependence of the produced gluon spectrum shows N paxt scal- 
ing in some limiting case, it can differ significantly from a simple N paTt or N co \\ scaling in 
the experimentally accessible regime depending on the shape of the gluon distribution . 
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A detailed comparison of perturbative saturation models to data not only requires 
the knowledge of the distribution and the improvement of the calculation of the gluon 
production yield beyond the factorized expression (J2J) used in this paper. It also requires 
the inclusion of fragmentation functions for gluons into pions, and a discussion of their 
possible medium-dependence [3]. 

On the qualitative level however we observe that the gluon distributions which lead 
to the Cronin effect in d-Au collisions also lead to the Cronin enhancement in the Au-Au 
collisions. And vice versa, if no Cronin effect appears in Au-Au, none is seen in d-Au 
collisions. Given the recent experimental observation of the Cronin enhancement in d-Au 
collisions at RHIC JH] this supports the view that significant final state ("quenching") 
effects are needed in order to account for the Au-Au data [TH] . 

Note added. When preparing the revised version of this paper we were made aware 
of [2D] and which also study the effects of saturation on the Cronin enhancement. 
These references agree with our results regarding the MV gluon. Regarding the evolved 
gluon, the detailed numerical study is reported in [TTj . 
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